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Abstract 

o : 

We present a systematic analysis for the rare baryonic exclusive decays of A& — > A/ + / (/ = 
Ph! e, [A, r). We study the differential decay rates and the di-lepton forward-backward, lepton 
polarization and various CP asymmetries with a new simple set of form factors inspired 
by the heavy quark effective theory. We show that most of the observables are insensitive 
to the non-perturbative QCD effects. To illustrate the effect of new physics, we discuss 
our results in an explicit supersymmetric extension of the standard model, which contains 
new CP violating phases and therefore induces sizable CP violating asymmetries. 



1 Introduction 



A priority in current particle physics research is to determine the parameters of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements |]] in the standard model (SM). 
Due to the CLEO measurement of the radiative b — > 57 decay 0, some interest has 
been focused on the rare decays related to b — ■> sl + l~ induced by the flavor changing 
neutral currents (FCNCs). In the SM, these rare decays occur at the loop level and 
depend on the CKM elements. In the literature, most of studies have been concentrated 
on the corresponding exclusive rare B-meson decays such as B — > K(*H + l~ [[|. However, 
these exclusive modes contain several unknown hadronic form factors, which cannot be 
measured in the present B-meson facilities unlike the kaon cases. Recently, we have 
examined the exclusive rare baryonic decays of A 6 — ► All (I — u, e, /i, r) j|, |5], [§] and 
found that some of physical quantities are insensitive to the hadronic uncertainties. 

In this paper, we give a systematic study on the baryonic decays of — > Al + l~ . We 
will explore various possible CP even and odd asymmetries to show how the hadronic 
unknown parameters are factored out in most of cases. To illustrate CP violating effect, 
we will also discuss an explicit CP violating model with SUSY. 

The paper is organized as follows. In Sec. 2, we give the effective Hamiltonian for 
the decays of A& — > All and the most general form factors in the A& — > A transition. In 
Sec. 3, we derive the general forms of the differential decay rates. In Sec. 4, we study the 
di-lepton forward-backward, lepton polarization and various CP violating asymmetries. 
We perform our numerical analysis in Sec. 5. We present our conclusions in Sec. 6. 



2 Effective Hamiltonian and form factors 

In the SM, the effective Hamiltonian for b — ► sl + l~ is given by 

(~i 10 

n = -A^v tb v: s Y,cMOM (i) 

where the expressions of the renormalized Wilson coefficients Cj(//) and operators Oi(n) 
can be found in Ref. J7|. From Eq. ([!]), the free quark decay amplitude is written as 



M{b^sin-) = ^v tb v* 



s [C e 9 ff 0) 1 ,P L - ^C< 7 e// 0) ia^PjA b T-fl 



V2: 

+sC 10 i^P L b h^l\ (2) 

with Pl(r) = (1 =F 75)/2. We note that in Eq. @, only the term associated with the 
Wilson coefficient C w is independent of the /1 scale. Besides the short- distance (SD) 
contributions, the long-distance (LD) ones such as that from the cc resonant states of 
...etc are also important for the decay rate. It is known that for the LD effects 
in the B-meson decays f| ||, 0, [11], [12], [R|, both the factorization assumption (FA) 
and the vector meson dominance (VMD) approximation have been used. In baryonic 
decays, we assume that the parametrization of LD contributions is the same as that in 
the B-meson decays. Hence, we may include the resonant effect (RE) by absorbing it 
to the corresponding Wilson coefficient. In this paper as a more complete analysis we 
also include the LD contributions to the decay of b — > S7, induced by the nonfactorizable 
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effects [H , [T5[ ]. The effective Wilson coefficients of C 9 
standard form 



eff 



and Cj can be expressed as the 
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Here s') combines the one-loop matrix elements and the LD contributions of operators 
Ox - Oq, C' b ^ sl is the absorptive part of b — > 57 [|T6j with neglecting the small contribution 

US ' 



m c /m b , s' = q 2 /ml, r] = a s (m w ) / ac s (/J,) , x t = mf/m^, Mj (Tj) 



arc 



from V ub V* s , z 

the masses (widths) of intermediate states, |a>| < 0.15 describing the nonfactorizable 
contributions to b — > 57 decay at q 2 = [13], |T5| , and the factors kj are phenomenological 
parameters for compensating the approximations of the FA and VMD and reproducing 
the correct branching ratios of B(A b -> AJ/^f -> A/ + /~) = £(A 6 — > AJ/\I/) x B(J/ty 
Z + /~) when we study the A& decays. We note that by taking ~ — l/(3Ci + C2) and 
-B(A{, — > AJ/ty) = (4.7 ± 2.8) x 10~ 4 , the kj factors in the A& case are almost the same as 
that in the B-meson one [Q. The Wilson coefficients (WCs) at the scale of \x ~ ~ 4.8 
GeV are shown in Table [l]. 



Table 1: Wilson coefficients for m 4 = 170 GeV, fi = 4.8 GeV. 







c 2 


c 3 


c 4 


c 5 


-0.226 


1.096 


0.01 


-0.024 


0.007 




c 6 




c 8 


c 9 




-0.028 


-0.305 


-0.15 


4.186 


-4.559 
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Using the form factors given in Appendix A, we write the amplitude of A b — > Al + l 

as 

M (A 6 - AZ+Z") = ^ + HSLl} (7) 



where 



#f = A^(A 1 P R + B 1 P L )A b + Aia^q l/ (A 2 P R + B 2 P L )A b , 
H% = A-f{D 1 P R + E 1 P L )A b + Ai(T^q u (D 2 P R + E 2 P L )A h (8) 
+q»A (D 3 P R + E 3 P L )A b , 

Ll = hsrut (9) 



with 



a _ n eff ft - 9i 2m b n eff fl + 9 T 

- °9 — -r L i ~Y~ 



fi+jh _ 2m brieff fT -gf 

2 q 2 



r> _ r<eff±i ill » n e 

°i — °9 2 

fi-9 



A = C 
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^ = C 10 ^^. (10) 

and i — 1, 2, 3. 

The processes for the heavy to light baryonic decays such as those with A5 — > A have 
been studied based on the heavy quark effective theory (HQET) in Ref. [|17[] and it is 
found that 

(A{p A )\sTb\A b (p Ab )) = u A (F 1 (q 2 )+ j6F 2 (q 2 ))Tu Ab (11) 

where T denotes the Dirac matrix, v = p Ab /M Ab is the four-velocity of A b , q = p Ab — p A 
is the momentum transfer, and F 12 are the form factors. Clearly, there are only two 
independent form factors F 12 in the HQET. Comparing with the general forms of the 
form factors in Appendix A, we get the relations among the form factors as follows: 

9i = fi = fl = gl = F x + VrF 2 , 

n r V rV F 2 

g2 = f 2 = g-i = fs = g T = Jt ~ 



M Ah 



g b T = /r = 0, 

T rT F 2 2 

^ = -^jh(MA b + M A ), /J — tt— (M Ab — M A ) , (12) 
M A b M Ab 

where r = M\jM\ h . From the CLEO result of R = -0.25±0.14±0.08 @, we know that 
1-^2 1 < l-^il- Due to Eq. (0), only fi (g x ) and /J (gj) are proportional to F x and therefore, 
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they are large, whereas all the others are small since they are related to the small form 
factor F2. Furthermore, from Eq. (fLO), we find that {/ T } and {g T } are associated with 
C7 which is about one order of the magnitude smaller than Cg and C\q so that their effects 
to the deviation of the results in the HQET are small. Hence, with the information of 
the HQET, we can make a good approximation for the general form factors of transition 
matrix elements given in Eqs. (|7|) and ([IT]). Altogether, we have the following relations: 



/ 

h -9i 
fi + 9i 

h + 9i <f 



h + 91 

2 ' fi + gi 

T T 

x 92 91 92 ■> 

"j r rT fT ' 

J2 Jl J2 

/2 + 92 

h + 91' 



(13) 



In the HQET, it is easy to show that 



Ma 



h + 92 
Ji+91 



Fi + ^F 2 ' 



(14) 



3 Differential decays rates 



In this section we first present the formulas by including the lepton mass for the double 
differential decay rates with respect to the angle of the lepton and the invariant mass of 
the di-lepton. In the following we only show the results of the SM with the form factors 
in Eq. ([131) . The general ones with including right-handed couplings are presented in 
Appendix B. 

Introducing dimensionless variables of t — p^ b ■ p\/M^ , r = M\/M\ . rhi = mi/M\ b , 
rhb = nib/M\ b , and s = q 2 /Ml b , the double partial differential decay rates for At, — > A Z + l~ 
(I = e, fj,, r) can be written as 



d 2 F G 2 F a 2 em \l 5 

dsdz ~ 768tt 5 Ab 



1 - 
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tRn (s,z) 



(15) 



where 



R A , (s, z) = I (s, z) + zh (s, z) + z I 2 (s, z) 



(16) 



and 



J (s, z) 
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+6m 2 t 
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(2m b pf C 7 e// 



(2m 6 p) 2 C e 7 ff 



+ 
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+ 



IC 
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10 1 
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m. 



+6Vr(l-i) 4 1 + 2-i 



+ps 
+ 12 I 1 
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,2\ 



U 9 



m fe p 



c 7 e// 



+ 1 



2— J- ] |Cio| 2 



2^- ] m 6 (t - r) (l + sp 2 ) HeCf Cf 



-12 ^1 + 2^f j m bv ^spReCfC 



/■toil* 



-6 



x 



s(l-t) (t-r)--((l-r) 2 - S 2 ) 



4m? 



C 7 e// " + sp 2 



-6m 2 (2r - (l + r)t) 



2m;, 



+ |C 10 | 

2 



\c 



10 



(17) 



MM) = 3^1-^0( S ){ S [l-2v^p-(l-r)p 2 ]ReCfC* lo 



+2m fe (l - sp 2 ) ReCfC* } , 



h (s, z) 



(2m 6 p) 2 |C 7 e// 



+ 



c Q e// 



|C 



10 



4m? 



C 7 e// " + sp l 



r eff 



+ IC; 
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(19) 



with z = pB ■ pi+ being the angle between the momenta of A& and l + in the di-lepton 
invariant mass frame and <p (s) = (1 — r) 2 — 2s (1 + r) + s 2 . Here, for the simplicity, we 
have not displayed the dependence of the p scale in effective Wilson coefficients. We note 
that the main nonperturbative QCD effect from / has been factored out in Eq. flT5]) . The 
function R\ b (s, z) is only related to the two parameters of 5 and p which become one in 
the HQET. Since p is the ratio of form factors and insensitive to the QCD models, the 
QCD effects in the baryonic di-lepton decays are clearly less significant. Therefore, these 
decay modes are good physical observable to test the SM. 

After integrating the angular dependence, the invariant mass distributions as function 
of s are given by 



ds 



384vr 5 Al 



d r(A^AW-)_ G | aL A? ^ 



\ 



1 - ^-p r Ai w 



where 



with 



R Ab (s) = r\ ( s ) + r 2 ( s ) + r 3 (s) 



(20) 



(21) 



6v^s -2m 6 p M + 2^J ReCf Cf * 
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off 
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(22) 



(23) 



(24) 



The limits for s are given by 



4m 2 < s < (l - v 7 ^) 2 • 



(25) 



From Eqs. (p2|) -(p4|), we see that p appears either as y/rp or p 2 which is small since 
r ~ 0.04 and |p| ~ 0.25. 



4 Lepton and CP asymmetries 
4.1 Forward-backward asymmetries 

The differential and normalized forward-backward asymmetries (FBAs) for the decays of 
A b -> Al+l- function of s are defined by 



dA 



FB (S 



ds 



1 , A rf 2 r(s,5) /-o 7 ^ rf 2 r( s ,z) 

— : — — / dz 



dsdz 



dsdz 



(26) 



and 



A FB (s) 



dT (s) /ds 



l d , d 2 T (s,z) 
o (isrfz 



l dsdz 



(27) 
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respectively. Explicitly, using Eq. flI5[), we obtain 



dA FB (s) _ G>1A* / 4 m[\ f ^ 



r/.s 2 8 tt 5 ^ '- u 

and 



/ \ 3 / , , . / 4m? Rfb (s) , \ 



where 



i2 F B(s) = s [l - - (1 - r) p 2 ] ReCf q o 

+2m b (l - sp 2 ) ReCf C* 10 . (30) 

It is known that the FBA is a parity-odd but CP-even observable, which depends on 
the chirality of the leptonic and hadronic currents. In order to obtain one power of z 
dependence, the related differential decay rate should be associated with TrL^L b v . This 
explains why the FBAs depend on ReCg^C^ and ReC^^ C^ . However, unlike that in the 
decays of B — > Kl + l~ where the FBAs are always zero since they only involve vector and 
tensor types of currents, the transition matrix elements in the baryonic decays preserve 
the chirality of free quark interaction. 



Similar to the B-meson decays [Q, pl| the FBA in Eq. (|29|) vanishes at sq which satisfies 
with the relation 

KeLg LiQ _ _________ Ke ^ C 10 . (61) 

We will see later that the vanishing point is only sensitive to the effects of weak interaction. 
4.2 Lepton polarization asymmetries 

To display the spin effects of the lepton, we choose the four-spin vector of Z + in terms of 
a unit vector, £, along the spin of l + in its rest frame, as 

s% = ^4 s + = i + W ^ r -P + , (32) 

mi E i+ + mi 

and the unit vectors along the longitudinal and transverse components of the l + polariza- 
tion to be 



P+ 

\p+y 

Pa x p 4 



\Pa x P+\ ' 

£n = e L x e T , (33) 



respectively. 



7 



Defining the longitudinal and transverse Z + polarization asymmetries by 



dT fe -i=l) -dT (ii ■ 1 = -1 
Pi (s) = ) % ( ^4 f , (34) 



with i = L and T, we find that 



\ 



Amf R L (s) 



f RaM 



-Trnii 



\ 



4mf 



- Rt(s) 



(35) 
(36) 



where 



Ri 



-ReCg ff C^ 



Ri 



1 - r) + s (1 + r) - 2s 2 + Qy^rps (1 - r + s) 
+p 2 s (2 (1 - rf - s (1 + r) - s 2 )" 
-6m 6 ReCf C* [(1 - r - s) (l + p 2 s) + 4v^psj , 
1 - 2^rp - p 2 (1 - r) 



ImCfCt 



+ ^(l-p 2 s) ImCfQ . 



(37) 



(38) 



Here we do not discuss the normal polarization (Pn) because the nonperturbative effects 
from the form factors are large at the small s region and moreover, the dependence of 
Wilson coefficients is similar to the invariant mass distribution 0. We note that the 
longitudinal lepton polarization of Pl in Eq. (|35"D is also a parity-odd and CP-even 
observable just like the FBA, whereas Pt in Eq. (|36|) a T-odd one which is related to the 
triple correlation of s + ■ (p\ x p + ). In general, Pt can be induced without CP violation as 
the cases in B-meson Et 



and kaon |21| decays. However, we expect that they are small. 
Moreover, such effects can be extracted away while we consider the difference between 
the particle and anti-particle as discussed in the next section. 



4.3 CP asymmetries 

In this subsection, we define the following interesting direct CP asymmetries (CPAs) by 

dT-dT 



A r 

A FB 



dT + dT ' 
dYpB — dYpB 
dT + dT 



dT (£■ 6i) - dT (£■ 6i) 

A - - { dT + dT > - L ' T 



(39) 
(40) 

(41) 



where we have used dT + dT as the normalization. The above four CPAs are CP-odd 
quantities and they are CP violating observables. For Ar t pB,p L in Eqs. d39|)-(^ID, to 
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display the difference of the physical observable between the particle and anti-particle, 
it is necessary to have the strong and weak phases simultaneously in the processes. In 
the decays of b — > sl + l~ (I = e, /i,t), the strong phases are generated by the absorptive 
parts of one-loop matrix elements in operators Oi ~ Oq and LD contributions. However, 
since Pt is a T-odd observable and only related to the imaginary couplings, even without 
strong phases, we still can have nonzero values of Ar(Aft — > A/ + /~). For b — > sl + l~ and 
b — > sl~l + decays the Wilson coefficients C 9 ^(p) and C 7 (pj) in Eqs. (H) and (f|) can be 
rewritten as 

C° 9 *(p)+iC^(ii), 
C° 7 (p) + iC? s (p), 

C° 7 *(p) +i C? s (p), (42) 



C 9 (p) + ReY(z,s') , 
C 7 (p,)+ReC' 7 (p,,q 2 ), 
ImY (z, s') , 

lmC' 7 (n,q 2 ), (43) 

where we have assumed that the strong phases are all from the SM and there are no weak 
phases in absorptive parts. We note that there is no strong phase in C\q. 

According to Eqs. (^), (p9|), (|35D, and (|36|), the CP asymmetries are all related to 
the following combinations: 



C e 9 ff (») = 

cv e// M = 
c e 7 ff (») = 

with 

c°M = 
c°M = 

ct(») = 

Cf s (p) = 



ReC e 9 ff C e 7 ff * - ReC e 9 ff C e 7 ff * 



ImC e 7 } JC* w -Im e 7 s iC- 



eff, 



°7,9 



°7,9 



2Cf s lmC^ + 2C? s ImC$ , 
2C!$ImC w , 

2/mC 7 ° 9 C , * + 2C^ b 9 s ImC w 



Explicitly, the CP asymmetries in Eqs. fl3D|), (plf), and ([IT]) are found to be 



(44) 



Ar 
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mt 



2^ps + {t - r) 1 + sp 2 ) C™ s ImC° 7 + Cf s ImC 



x 



+ 



X 



X 



/ mt 

,2 „2siabs 



At 2 (l-^r) + 3(l + r)t + Qmft 



4m'ip z C^ s lmC° 7 + C^ImCg 



~tabsn 



-2t 2 1 + 2- 



mt 



Ami 



- 4r 1 - 



Gf s ImCl^ + sp 2 Cf s ImC 



m. 



mt 



+ 3 (1 + r) t - 6— (2r - 1 - tr) 

s 



mt 



4mgCf s /mC 7 u + sC^ImCl 



abs : 



&\frs<5 1 



mf 



Cf s ImCl 



(45) 
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A 
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1 - 



4m? 



2RA b (s)\ q 
+ 2m b (l - sp 2 ) C« 

1 _ Am} 



s)ImC 10 [s (l - 2^p - (1 - r) p 2 ) C, 



abs 
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abs 



+s (l - sp 2 ) (1 + r) - s 2 
+Qm b Cf s [(1 - r - s) (l + sp 2 ) + A^ps] } , 

3%ifii . Am 2 



4#A 6 (*)\ q 



1 - 



s0 (s) { (lmC°C* + Cf s ImCia 



x (l - 2v^p - (s + 2t - 2r) p ; 



sp 



[lmJC^C{ Q 



Cf s ImC w 



(46) 



1 - -j±IrnC 1Q {C« bs [6V^ps (1 - r + s) + (l + 2sp 2 ) (1 
2 + sp 2 ) 



(47) 



(48 



As seen from the above equations, Ap is related to IH1C7 and ImCg, while App, Ap L and 
Ap T depend on ImC\Q. Moreover, for small values of Cg bs ImCio and C % bs I mC 10, Ap T 
would still be sizable because ImCgCl or ImC®C* would be large. 



5 Numerical analysis 

In our numerical calculations, the Wilson coefficients are evaluated at the scale p ~ m b 
and the other parameters are listed in Table 1 of Ref. 0]. From Eq. (13), we know that 
the main effects to the deviation of the HQET are from 5. By using a proper nonzero value 
of 5, we will see later that the deviations of the decay branching ratios of A b — > hl + l~ 
are only a few percent. Since there is no complete calculation for the form factors of the 
A b — > A transition in the literature, we use the form factors derived from QCD sum rule 
under the assumption of the HQET, given by 

fm 2 ) = nr^xTTi > ( 49 ) 

1 + aq z + og 4 

with the parameters shown in Table 1 of Ref. || . In order to illustrate the contributions 
of new physics, we adopt the results of the generic supersymmetric extension of the SM 



22|j in which 

C* USY = -1.75 {5^) LL - 0.25 (6% 3 ) LR - 10.3 

23) LR ' 



Cl USY = 0.82 {51 



Cf USY = -9.37 (8&) LR + 1.4 (6^) LR (6 U 33 ) RL + 2.7 (S^) LL , (50) 
and take the following values instead of scanning the whole allowed parameter space: 

(*sb)m ~ 0.65, 



W M - 3X10-V-, 
Wlr ~ -O-Se'S, (51) 
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where (<5?)ab (hj = 1, 2, 3 and A, B = L, R) denote the parameters in the mass insertion 
method, which describe the effects of the flavour violation. The set of the parameters in 
Eq. (|5lD satisfies with the constraint from B — > X s 7 on C 7 = C* M + C* USY ||. Hence, 



the numerical values of the SUSY contributions to the relevant Wilson coefficients are as 
follows: 

ReC^ USY ~ 0.06, ImC^ USY ~ -0.29, 
ReCl USY ~ -0.46, ImCl USY ~ -0.46, 
ReCff 3 * ~ 4.78, JmCf c/5y ~ 4.50. (52) 

We note that the contributions of the minimal supersymmetric standard model (MSSM) 
to 6 — > sZ + /~ can be found in Refs. and |2~| . 

To show the typical values of various asymmetries, we define the integrated quantities 

_ J t> max 

Q = Q(s)ds (53) 

J s min 

where Q denote the physical observables with s min = 4m; and s max = (1 — \Zr) 2 . 



5.1 Decay rates and invariant mass distributions 

We now discuss the influences of 5, p, and ui on the branching ratios (BRs) of — > Al + l~ 
decays in detail. The effects of kj for compensating the assumption of the FA and VMD 
have been analyzed in 0. In Table |2|, we show the BRs by choosing different sets of 
parameters. Our results are given as follows: 



Table 2: BRs (in the unit of 10 6 ) for various parameters with uj = and neglecting LD 
effects. 



Parameter 


A b -> Ae+e" 


A& — ► Ap + p 


A b —> At + t~ 


HQET 


2.23 


2.08 


1.79 x lO" 1 


5 = 0.05 


2.36 


2.21 


1.86 x 10" 1 


5 = -0.05 


2.09 


1.96 


1.71 x lO" 1 


p = 0,5 = 


2.52 


2.38 


2.66 x lO" 1 


C 7 = 0,5 = 


2.36 


2.34 


2.23 x lO^ 1 


C 7 = -C$ M ,8 = 


3.34 


3.19 


2.76 x 10- 1 



1. By taking \5\ = 0.05 which means 10% away from that in to the HQET, we clearly 
see that the deviations of the BRs are only 4 — 6%. It is a good approximation to 
neglect the explicit 5 term in Eqs. fllTD , (f2~2|) and (flop. Hence, / = (fx + gi) /2, 
which also owns the 5 effect, is the main nonperturbative part. 

2. If p = 0, the effects are about 10% for e and p modes but 48% for r one. 

3. If one neglects the contribution from C 7 , the influences on B(A b — > Al + l~) for e, 
p and t modes are about 5%, 12% and 24%, respectively. However, taking the 
magnitude of C 7 is the same as the SM but with an opposite sign, the deviations 
are all over 50%. 
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The contributions of the parameter uj to the BRs of A& — > Al + l~ are listed in Table 
|3] and the invariant mass distributions are shown in Figure 1. From Table |3], it is clear 
that the nonfactorizable effects are small on the BRs. However, those directly related to 
uj effects such as Pt and CP asymmetries will have large influences. 

Table 3: BRs (in the unit of 10~ 6 ) without LD effects with different values of uj 



Mode 


A b -> Ae + e" 


Ab — > A/i + /i 


At -> Ar+r~ 


u = 0.15 


2.24 


2.12 


1.89 x 10" 1 


u = 0. 


2.23 


2.08 


1.79 x 10- 1 


uj = -0.15 


2.25 


2.06 


1.71 x 10" 1 



As for the new physics contributions, using the values of the SUSY model in Eq. (|5"i2"D, 
we show the results in Table f|. Although the deviations of the BRs to the SM are not 
significant, they have a large effect on the lepton and CP asymmetries which will be shown 
next. 



Table 4: BRs (in unit of 10 6 ) in the generic SUSY model. 



Model 


A;, -> Ae+e~ 


Afe — »■ Ap + p 


Afe -> Ar+r~ 


SUSY 


2.47 


2.24 


1.79 x 10- 1 



5.2 Forward-backward and lepton polarization asymmetries 



From Eq. and R = F 2 j F x ~ -0.25 in the HQET, we have that p ~ -0.26. We note 
that p is defined by the ratio of the form factors and it is expected to be insensitive to 
the QCD models. With s max ~ 0.64, we obtain s max p 2 ~ 0.04, (1 - r)p 2 ~ 0.06 and 2^/rp 
~ 0.2. Using these values, one can simplify Eqs. ( PU[ ) and Q3"T| ) to 



R 



FB (S) 



s (l - 2v^p) ReCf C* + 2m 6 ReCf C* 10 



and 



ReCf C* ~ 



2rhh 



(54) 



(55) 



s (l-2VFp) 

respectively. It is easy to see that so is only sensitive to the Wilson coefficients. The result 
is similar to the case in B — > K*l + l~ decay || where the approximation of the large 
energy effective theory (LEET) [^3J is used. As for the lepton polarization asymmetries, 
with the same approximation, Eqs. fl37D and (B3) can also be reduced to 



Ri 



Ri 



-ReCf Cio [l + s - 2s 2 + 6v^ps (1 + s 
-6m 6 ReCf C* 10 [l - s + Ay/rpa 

1 - 2v^p) ImCf C* + ^ImCf C* 



10 ' 



(56) 
(57) 
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respectively. Hence, the lepton asymmetries are all more sensitive to the Wilson coeffi- 
cients than the nonperturbative QCD effects. 

It is worth to mention that the effects of u, introduced for the LD contributions to 
b — > 57 and absorbed to , will change ReC| ff in the SM such that So is also shifted. 
Therefore, in terms of Sq, we can also theoretically determine to by comparing the result 
with that of u = 0. Another interesting quantity is T-odd observable of Pt which is 
proportional to CioImC^ in the SM. Due to the enhancement of Cxo, a nonzero value of 
10 will modify Pt enormously. As for the other asymmetries, the effects are insignificant. 
The estimations of integrated lepton asymmetries with different values of u in the SM are 
displayed in Table ||] and the corresponding distributions are shown in Figures 2 — 4. 

Table 5: Integrated lepton asymmetries in the SM without LD effects. 



Parameter 


Mode 


10 2 A FB 


W 2 Pl 


10 2 P T 


u = 0.15 


Ab — > A/x + /x~ 


-14.37 


59.50 


0.11 


A b -> Ar + r~ 


-3.98 


10.70 


0.53 


u = 0. 


A b -> A/X+/X" 


-13.38 


58.30 


0.07 


A b -> At + t~ 


-3.99 


10.84 


0.39 


u = -0.15 


Ab — > An + n~ 


-12.24 


56.70 


0.04 


A b -> At + t~ 


-4.00 


10.94 


0.23 



To illustrate the new physics effects, the integrated lepton asymmetries in the generic 
SUSY model with u> are listed in Table || and their distributions as a function of q 2 /Mx b 
are shown in Figures 5 — 7. From the figures, we see that SUSY effects make the shapes of 
lepton asymmetries quite differ from that in the SM. We summary the results as follows: 

1. Since the SD contributions to CgC% and i?eC 7 C* are —1.40 and —1.35, respectively, 
which violate the condition in Eq. fl55]) , the vanishing point is removed. 

2. Due to the factor of rh b /s, from Figure 7, we see that ImC^ C* Q has a large effect 
on Pt in the small s region. 

3. In the SUSY model, Pt could reach 1% and 10% for the light lepton and r modes, 
which are only 0.2% and 3% at most in the SM, respectively. 



Table 6: Integrated lepton asymmetries in the generic SUSY model with u = 0. 



Mode 


10 2 A FB 


10 2 P L 


10 2 P T 


A b — > Afi + fi~ 


-10.53 


24.46 


-0.57 


A b -> At+t" 


-1.84 


4.40 


-2.51 
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Table 7: CP asymmetries in the generic SUSY model for different values of uj 



Parameter 


Mode 


10 2 A r 


10 2 A FBA 


10 2 A Pi 


10 2 A Pt 


uj = 0.15 


A5 — > A/i + /i~ 


2.05 


-2.62 


6.48 


-0.53 


A 6 -> Ar+r~ 


1.83 


-0.79 


1.94 


-2.01 


uj = 0. 


A5 — > A/i + /i~ 


1.59 


-1.89 


5.00 


-0.47 


A b -> Ar+r~ 


1.38 


-0.59 


1.53 


-2.21 


uj = -0.15 


At — > An + n~ 


1.05 


-1.06 


3.34 


-0.40 


A b -> Ar+r" 


0.89 


-0.37 


1.06 


-2.41 



5.3 CP asymmetries 

In the SM, for b —>■ sl + l~ , the relevant CKM matrix element is V t bV t * which is real under 
the Wolfenstein's parametrization. Nonzero CPAs will indicate clearly the existence of 
new physics. We remark that the CPAs can be in fact induced by the complex CKM 
matrix element V u bV* s which is also the source of the direct CPA in B — > X s 7 in the SM. 
However, we expect that such effects to the CPAs in b — > sl + l~ are smaller than that in 
B — > Xg'-f where the CPA is less than 1%. The main reason for the smallness is because 
of the presences of Cg and C\q contributions to the rates of b — > sl + l~ , which are absent 
in B -> X,7. 

With the values in Eq. ([52]), the averaged CPAs in the generic SUSY model for 
Aft — > A/ + /~ are listed in Table |7] and their distributions as a function of s — q 2 jM\ are 
shown in Figures 8 — 11. The results are given as follows: 

1. From Eqs. (flop and ([48]) , we see that the terms corresponding to C^ImC® and 

IuiCjCIq + C^ bs ImCio are associated with a factor of rhb/s. If sizable imaginary 

parts exist, in the small s region the distributions will be significant. Due to this 

reason, in Figure 8a one finds that Ar(s) for A& — > Al + l~ (/ = e, /x) increase as s 

decreases. On the other hand, if the term with rhb/s in Eq. (f48|) is dropped, the 

distributions of Ap T (s) for e and /1 modes do not contain zero value. We note that 

with the values in Eq. ([52"D, the main effect on Ap(s) in the small s region is from 
ni 

2. Ap L (s) for all lepton channels and Ap T (s) for the r one could be over 10%, while 
the remaining CP asymmetries are at the level of a few percent. We remark that if 
we can scan all the allowed SUSY parameters, the asymmetries except Ap T (s) for 
lighter lepton modes would reach up 10%. 

3. It is known that Ap T (s) is a T-odd observable and the other CPAs belong to the 
direct CP violation which needs absorptive parts in the processes. This is the reason 
why the distributions of Ar(s), Afba{s) and Ap L (s) around the RE region have 
the similar shapes but are different from that of Ap T (s). Moreover, all the direct 
CPAs are sensitive to uj unlike the cases of the CP conserving lepton asymmetries 
discussed in Sec. 5.2. 
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6 Conclusions 



We have given a systematic study on the rare baryonic decays of A& — > A/ + /~ (/ = e, p, r). 
For the — > A transition, we have related all the form factors with F\ and F2, and we 
have found that 5 = and p ~ R = F2/F1, in the limit of the HQET. Inspired by the 
HQEF, we have presented the differential decay rates and the di-lepton forward-backward, 
lepton polarization and four possible CP violating asymmetries in terms of the parameters 
/, 5 and p. We have shown that the non-factorizable effects for the BRs and CP-even 
lepton asymmetries are small but large for and the direct CPAs. We have also demon- 
strated that most of the observables such as Afb, Pl,t and A a (a = T, FB, P L andP^)- 
are insensitive to the non-perturbative QCD effects. We have illustrated our results in 
the specific CP violating SUSY model. We have found that all the direct CP violating 
asymmetries are in the level of 1 — 10%. To measure these asymmetries at the no level, 
for example, in the tau mode, at least 0.5n 2 x (10 9 — 10 10 ) A& decays are required. It 
could be done in the second generation B-physics experiments, such as LHCb, ATLAS, 
and CMS at the LHC, and BTeV at the Tevatron, which produce ~ 10 12 66 pairs per year 
26| Finally we remark that measuring these CPAs at a level of 10~ 2 is a clear indication 



of new CP violation mechanism beyond the SM. 
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Appendix 

A. Form factors and decay amplitudes 

For the exclusive decays involving A b (p Ab ) — > A(p A ), the transition form factors can 
be parametrized generally as follows: 



(A\s 7m b\A b ) 
(A\s 7m 75 6 | At) 
(A | sio^b \A b ) 

(A\ sia^^b \A b ) 



fiUKl^u^ + f 2 u A ia^ q v u Ab + foq^u^u^ (58) 

9iUAlvl5U Ab + g 2 u A ia lxv q v ^u Ab + g^u^u^, (59) 

fru A ia^u Ab + fr u A (^q„ - ^ u q^) u Ab 

+/t ( P hQ» - p y%) u A u Ab , (60) 

g T u A ia^ 5 u Ab + g?u A (7^ - 7^) ^ 5 u Ab 

+9t { p n<lv ~ PvQn) u A ^ 5 u Ab , (61) 



where P = p Ab + p A , q = p Ab — p A and form factors, {/$} and {gi}, are all functions of 
q 2 . Using the equations of the motion, we have 



(M A + Ma J u A -f^u Ab 
(M A - M A J u A <y^ 5 u Ab 



(PA b + Pa)^ u A u Ab + iu A cr^ u q u u Ab , 
(?A b + Pa)^ u A -f 5 u Ab + iu A a^q v ^u Ab 



The form factors for dipole operators are derived as 



(A| sia^ v q v b \A b ) 
(A| sia^q^^b \A b ) 



fi' u AlfM u A b + fTu A ia^q"u Ab + f!q^u A u Ab , 
glu A ^^u Ab + #|Wo><f 75«A 6 + glq^u A ^u Ab 



(62) 
(63) 



(64) 
(65) 



with 



fT = [f¥ + f$(M A + M Ab ] 



fT 



{M Ab -M A y- 



S„2 



9 2 = 9t-9tQ' 

9t + 9t {M a — M Ab ) 



T 

9i 



T 

9i 



q 



(M Ab + M A ) 



T 

9 3 ■ 



(66) 



We now give the most general formulas by including the right-handed coupling in 
the effective Hamiltonian with a complete set of form factors. The free quark decay 
amplitudes for b — > sl + l~ are given by 



H (b -> sl + l~) 



Gfo:, 



em \r t/* 

v^r tbts 



(cip L + c*p R ) b h,i 



(c^P L + C*P R ) b 



2m 



-sia^ (C 7 L P R + C 7 R P L ) b 



(67) 
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where C[ and Cf (i = 7, 9, 10) denote the effective Wilson coefficients of left- and right- 
handed couplings, respectively. With the most general form factors in Eqs. (|58|) , (|59|). 
(H) and (H), and 
for the decays of 

M (A 6 -> A/ + / 

+Ai<vg" (A 2 P fl + P 2 P L ) A 6 Z 7m Z 



and (|65|) , and the effective Hamiltonian in Eq. (|67D, the transition matrix elements 
Al + l~ are expressed as 

= C ^^V tb V*{[A^(A 1 P R + B 1 P L )A b 



where 



+ 



A T/i (DiPfl + PxPl) A fe + Aia^g" (D^Pr + E 2 P L ) A b 



+q^K (D 3 P R + E 3 P L )A t 



(6* 



c, 



R 



fi + g t 2m b R fj - gf „ L fi~ 2m b L ff + g 



in 

,2 °7 



2 °7 



a 



L fi + gi 2m b L ff-gJ 



q 2 ° 7 2 



2 g 
fi-gi 2m, R fi+g. 



fi + 9i 



n R Ji ~T yi n L fi 9i 

°10 o r O 10 - 



n L fi + gj _,_ ^fl /i gi 

°10 7, r °10 o 



9 2 



(69) 



B. Differential decay rates 

Using the transition matrix elements in Eq. (|68|) , the double differential decay rates 
can be derived as 



where 



with 



dT 

dsdz 



H 2 rv 2 \ 2 



7687T 5 



4mf 



f 2 RA b (S,Z) 



R Ah (s, z) = I Q (s, z) + zli (s, z) + z 2 I 2 (s, z) 



(70) 



(71) 



In (s, z) 



—Q\/rs 



m 



m. 



1 + 2-f ReAiB* + 1 - 6-f ReDxE* 



q- 



r — s 



+6mft(|A 1 | 2 + |P 1 | 2 -|P 1 | 2 -|P 1 
+ 12m 2 M Ai)V ^ (1 - t) (ReDiDg* + ReEiE; 
+ 12m 2 M Ab (t - r) (ReDxE 3 * + ReD 3 E'*) 
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h 0, z) 



h (s, z) 



+6M AbV ^ S (l-t) 



1 + 2^- J (ReAiA^ + ReBxBg) 



+ 1- 2- 



mr 



q 2 J 



(ReDxD; + ReEiE^ 



-6M Ab s(t-r) 



1 + 2-f (ReAiB^ + ReA 2 B*) 



m 



+ 1 - 6-i (ReDiEg + ReD 2 E*) 



-6Mi^ 2 1 + 



2-1 J ReA 2 B^ - 6Mi bV ^s 2 ( 1 - 6-f J ReD 2 E; 



s (1 _ t ) (t - r) - - ((1 - rf - s 2 )\ (\A 2 \ 2 + \B 2 \ 2 + \D 2 \ 2 + \E 2 

-GM'im 2 (2r - (1 + r) t) (\A 2 \ 2 + \B 2 \ 2 - \D 2 \ 2 - \E 2 \ 2 ) 
+ 12m 2 Ml b st (ReD 2 D 3 + ReE 2 E 3 ) + 12m? M^y/rs (ReD 2 E 3 + ReD 3 E 2 ) 
3s0 (s) {- (ReAiD* - ReBiE*) + M Ab [^(ReA^ - ReBiE^) 

+ (ReAiE^ - ReBxD^) + ^ (ReA 2 D* - ReB 2 E*) - (ReA 2 E^ - ReB 2 D*) 
+M Ab (1 - r) (ReA 2 D; - ReB.E^)} , 



— A <P (s) [1 ~ 4^J [\Ai\* + + |£>i| + l^i 
+^Mi fc (a) (l - 4^ (|A 2 | 2 + |5 2 | 2 + |D 2 | 2 + \E 2 \ 



(72) 



where D' 3 — D 3 — D 2 and E' 3 = E 3 — E 2 , and z = p B ■ p t + denotes the angle between the 
momentum of A& and that of l + in the di-lepton invariant mass frame. 

C. Forward- backward and lepton asymmetries 



From Eq. (|S8"D, the functions of R Ab and Rfb m the differential and normalized FBAs 
for A 6 -> A/+Z" in Eqs. @ and (j29|) are given by 



1 r 1 



(73) 



and 

Rfb (s) 



-Re (AiDJ - BiE*) + M Ab [v^Re (AiD 2 - BiE 2 ) + Re (AxE 2 - BiDJ 
+v^Re (A 2 D* X - B 2 E*) - Re (A 2 E* - B 2 D*) 
+M 2 (l-r)Re(A 2 D* 2 -B 2 E;), 



(74) 

respectively. We can also define the longitudinal and transverse lepton polarization asym- 
metries. Explicitly, by the definition of Eq. (B3), we get 



Pt, 



I L 4:mf 



RaM 



1 l - { (s (1 + r - s) + (1 - r) 2 - s 2 ) Re (A X D* + B^) 
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-6s\/r 



-8(1 



2(1 



r? - 2s 2 



Re (A 2 D* 2 + B 2 E;) 



Re (AiE* + BiDJ) + M^sRe (A 2 E^ + B 2 D;; 
+3sM A (1 - r + s) [Re (A^* + BiE*) + Re (A 2 D* + B 2 E*)] 
-3sM Ab (1-r-s) [Re (AiE*. + B^) - + e (A 2 E* X + B 2 D*)] , 

3 , 



-nuii 



1 - 



4mf 



1 



{— Im (AiD* — BiE*) 



+M A [Im (AxD; - B 1 E* 2 ) + Im (A 2 D* - B 2 E*)] 
+M Ab [Im (A1E2 - BxDj) - Im (A 2 E* X - B 2 D*)] 
+Ml (1-r) Im(A 2 D;-B 2 E;). 
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Figure Captions 

Figure 1: BRs as a function of q 2 /M\ for (a) A b — > A/i + /i~ and (b) A 6 — »■ Ar + r~. 

The curves with and without resonant shapes represent including and no LD 
contributions, respectively. The dashed, solid and dash-dotted curves stand 
for uj = 0.15, 0, and —0.15, respectively. 

Figure 2: Same as Figure 1 but for the FBAs. 

Figure 3: Same as Figure 1 but for the longitudinal polarization asymmetries. 

Figure 4: Same as Figure 1 but for the transverse polarization asymmetries. 

Figure 5: FBAs in the generic SUSY model as a function of q 2 jM\ for (a) A5 — > A/i + /i~ 
and (b) A b — > At + t~. The solid and dashed curves stand for the SM and 
SUSY model, respectively. 
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Figure 1: BRs as a function of q 2 /M\ b for (a) A b — »■ Ajj, + jj,~ and (b) A b — »■ Ar + r~. The 
curves with and without resonant shapes represent including and no LD contributions, 
respectively. The dashed, solid and dash-dotted curves stand for u = 0.15, 0, and —0.15, 
respectively. 
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Figure 2: Same as Figure 1 but for the FBAs. 
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Figure 3: Same as Figure 1 but for the longitudinal polarization asymmetries. 
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Figure 4: Same as Figure 1 but for the transverse polarization asymmetries. 



26 





Figure 5: FBAs in the generic SUSY model as a function of q 2 /M? for (a) A5 — > Afi + fi~ 
and (b) A5 — > Ar + r _ . The solid and dashed curves stand for the SM and SUSY model, 
respectively. 
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Figure 6: Same as Figure 5 but for the longitudinal polarization asymmetries. 
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Figure 7: Same as Figure 5 but for the transverse polarization asymmetries. 
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Figure 9: Same as Figure 5 but for A FBj y. 
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